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ABSTRACT 


Honthly  Bean  values  o f  temperature  output  for  a  ten-year 
period  (1969-78)  from  a  three-dimensional  primitive  equation 
ocean  model  are  used  to  compute  isotherm  displacements  at 
eight  selected  grid  points  in  the  North  Pacific  Ocean.  The 
ocean  model  input  parameters  are  climatological  radiation, 
clouds,  surface  air  tempejsture,  surface  humidity  and  six- 
hourly  synoptic  winds  taken  from  Fleet  Numerical 
Oceanography  Center  analyses.  The  first  five  baroclinic 
modes  are  computed  from  the  mean  temperature  profiles  a.  -he 
eight  selected  model  grid  points  and  these  modes  are  then 
fit  in  the  time  domain  to  the  isctherm  displacements. 

The  fit  shews  evidence  of  baroclinic  free  Rossby  waves 
in  the  time  demain  at  the  selected  grid  points  south  of  40N 
and  indicates  the  predominance  of  the  first  baroclinic  meda. 
At  the  selected  grid  points  north  of  40N,  the  results  are 
inconclusive.  This  is  due  to  the  longer  periods  required  at 
the  northern  latitudes,  the  neglect  of  salinity,  and  mask_ng 
by  a  strong  seasonal  signal  caused  by  surface  haat^ng  and 
cooling  at  higher  latitudes. 
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INTRODUCTION 


Oceanic  general  circulation  models  have  been  under 
development  for  the  past  decade  and  a  half.  Evaluation  of 
these  newly  developed  models  was  initially  hindered  by  a 
lack  of  adequate  observational  data  from  ocean  regions. 
However,  ever  the  past  eight  years  high  quality  oceanic  data 
and  data  analysis  studies  from  portions  of  the  ocean  have 
been  made  available  to  the  large-scale  ocean  modeller.  Of 
particular  importance  and  interest  in  this  study  has  been 
the  analyses  of  temperature  data  resulting  from  the  "North 
Pacific  Experiffent"  (NORPAX)  begun  during  the  last  decade. 

Eased  upon  an  analysis  of  TRANSPAC  data  (a  sub-program 
of  NCRPAX  confined  to  a  region  3  0-50N  and  170E  -  150K),  Kang 
and  Bagaard  (1980)  determined  that  annual  temperature  fluc¬ 
tuations  in  the  area  3Q-40N,  16 OE  -  130W  consisted  of  first 
mode  Rossfcy  waves  with  random  phase  traveling  northwest  with 
a  wavelength  of  300  km.  For  other  latitudes  of  the  North 
Pacific  Ocean,  Prica  and  Saqaard  (1980)  determined  that 
first  mode  baroclinic  Rossby  waves  with  dominant  period  of 
seven  to  ten  years  existed  in  the  region  40-50N.  It  was 
also  determined  that  Rossby  waves  of  annual  period  did  not 
play  a  role  in  the  latitude  belt  40-50N  in  the  North  Pacific 


Ocean.  Magaard  (1933)  ,  ir.  a  paper  discussing  banclir.ic 
Rossty  wave  energetics,  summarizes  these  earlier  studies  by 
stating  that  first  mode  baroclinic  Rossby  waves  dominate  the 
internal  temperature  fluctuations  in  the  region  20-30N  ir. 
the  North  Pacific  Ocean  for  a  whole  range  of  periods  between 
five  months  and  ten  years.  3etween  30-40N  only  annual  or 
near-annual  first  mode  baroclinic  Rossby  waves  are  evident. 
It  is  the  results  of  these  studies  that  have  motivated  this 
investigation,  which  is  essentually  an  evaluation  of  whether 
a  general  circulation  ocean  model,  driven  by  observed  synop¬ 
tic  winds,  can  simulate  the  above  characteristics  of  the 
North  Pacific  Ocean  Rossby  wave  field.  Following  a  proce¬ 
dure  outlined  by  Emery  and  Magaard  (1976)  ,  this  study  exam¬ 


ines  the  possible  existence  of  baroclinic  Rossby  waves 
inferred  from  temperature  fields  produced  by  the  ocean 
model. 


II.  MODEL  DESCRIPTION  AND  ANALYSIS  PROCEDURE 

A.  MODEL  DESCRIPTION 

The  model  used  in  this  study  is  a  three-dimensional 
primitive  equation  model  developed  at  the  Naval  Postgraduate 
School  and  described  by  Haney  (1980) .  A  flat  bottom  and 
rigid  lid  are  prescribed.  The  model  is  driven  by  seasonally 
varying  (monthly  climatology)  solar  radiation,  clouds,  sur¬ 
face  air  temperature,  and  surface  relative  humidity.  The 
surface  winds  utilized  are  taken  from  six-hourly  synoptic 
analyses,  1969-1  978,  irade  at  the  Fleet  Numerical 
oceanography  Center,  Monterey,  California.  Salinity  is  neg¬ 
lected.  The  model  domain  shown  in  Figure  1  is  a  rectangular 
basin  extending  from  the  Equator  to  65N  and  from  145E 
12 5w .  The  horizontal  grid  spacing  is  approximately  two 
degrees  in  latitude  and  three  degrees  in  longitude.  The 
model  is  4000  m  deep  with  20  irregularly  spaced  levels,  14 
of  which  are  above  1000  m.  The  boundary  conditions  pre¬ 
scribed  at  the  bottom  acd  sides  are  rhose  of  an  enclosed 
thermally  insulated  basin. 
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in 


145  E  160  E  175  E  170  W  155  H  140  M  125 

Horizontal  domain  of  temperature  field  used  for  analysis  is 
shown  by  enclosed  region.  Triangles  mark  grid  points 
selected  fcr  Rcssby  wave  analysis. 

Figure  1.  Horizontal  Domain  of  Ocean  Model 


The  model  also  includes  a  dynamic  adjustment  mechanism 
as  described  by  Adamec,  et.  al.  (1981).  Essentially,  the 
mechanism  assures  that  if  the  local  gradient  Richardson  num¬ 
ber,  Ri,  falls  below  a  critical  value  (0.25)  then  the  temp¬ 
erature  and  the  velocity  components  at  the  two  adjacent 
levels  are  mixed  vertically  so  that  heat  and  momentum  are 


conserved,  tke  mixing  ratios  for  heat  and  momentum  are 
equal,  and  after  mixing  Ri  equals  the  critical  value. 

The  model  was  run  for  a  30- year  period,  with  the  first 
two  decades  considered  as  a  Mspin-up"  period  and  the  third 
decade,  1969-78,  considered  as  a  simulation  of  actual  condi¬ 
tions.  The  purpose  of  this  thesis  is  to  analyze  the  temper¬ 
ature  variability  simulated  by  the  model  during  this  last 
decade.  The  ccal  is  to  see  if  the  variability  is  consistent 
with  the  dispersion  relaticn  for  baroclinic  Rossby  waves. 

3.  ANALYSIS  EBOCEDORE 

1.  Governing  Equations  for  Baroclinic  Rossby  Waves 

The  governing  equations  for  baroclinic  Rossby  waves 
are  derived  from  the  frictionless  momentum  equations,  the 
hydrostatic  balance  and  the  continuity  equation  for  an 
incompressible  fluid  in  the  following  form: 


u.  +  ( V*  7)u  -  fv  +  ^  p  =  0  ,  (la) 

L  P  X 

vt  +  ( V*  V ) v  +  fu  +  i  py  =  0  ,  (lb) 

7»  V  =  0  ,  (lc) 

P_  +  Pg  =  0  »  (ld) 

Pt  +  (V«V)p  =  0  ,  (le) 
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where  V  is  the  velocity  vector,  u,  v  and  w  are  the  velocity 
components  in  the  eastward  (x) ,  northward  (y)  ,  and  upward 
(z)  direction,  respectively,  p  is  the  density,  p  is  the 
pressure,  g  is  the  acceleration  due  to  gravity,  f  is  the 
Coriolis  parameter,  t  is  time  and  7  =  )  •  The  sub¬ 

scripts  denote  partial  differentiation  with  respect  to  the 
subscripted  variable. 

The  vcrticity  equation  is  formed  using  equations 
(la)  and  (1b).  Subsequent  application  of  the  guasi- 
geostrophic  approximation  and  linearization  leads  to  the 
following  fcrm  of  the  vorticity  equation: 


(2) 


where  p'  and  w'  are  the  perturbation  pressure  and  vertical 
velocity  component,  0  =  df/dy,  p„  is  a  constant  reference 
density,  and  the  other  variables  are  as  defined  before. 
Additional  assumptions  used  in  the  above  linearization  pro¬ 
cess  are  that  the  lateral  shear  of  the  basic  flow  is  neg¬ 
lected  and  the  basic  state  horizontal  velocity  is  at  rest 
(i.e.,  u  *  v  *  0)  .  Prom  this  assumption,  equation  (1c)  and 
the  boundary  condition  that  w  =  0  at  z  =  0  (sea  surface)  it 
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follows  that  5=0.  Equation  (2)  contains  two  unknown  vari¬ 
ables,  p'  and  w*,  therefore,  an  additional  equation  is 
required. 

Substituting  from  (Id)  into  (1e)  after  multiplying 
by  ~g/P*  and  subsequent  linearization  yields  the  second 
equation, 

p'  +  p  N^w'  =0  ( 3a) 

zx  o 


or 


w’ 


=  -  ( 


2 

N  p 


(3b) 


where  N2  =  -(g/pJPz  is  the 
assumed  to  be  a  function  of 
using  (3b)  yields, 

(Pi„  *  Pjy)t  Bp'  *  f2!-^ 


Brun t-Vaisala  frequency  squared, 
z  only.  Eliminating  w'  in  (2) 


*zt'z 


=  0 


(4) 


Using  the  separation  of  variables  technique  the  perturbation 
pressure  is  written. 


p'(x,y,z,t)  =  Pn(z)  pn<x,y,t)  , 


(5) 
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where  n  is  the  acde  nuaber.  Substituting  (5)  into  equations 
(4)  and  <3fc)  gives. 


1 


d  .  1 

st7 


dP. 


n 


dz 


]  = 


?2(Pn>t  +  =  .  2 

r2,  *  v  n 


(6) 


and 


N7  ^ 


Pp  w' 
<Pn>t 


where  V2-  and  X*  is  a 

dinensions  cf  inverse  velocity. 

(6)  with  respect  to  z  yields. 


(7) 

separation  constant  with 
Differentiating  equation 


ri.  d?m 

V  ^ 


+ 


2  15, 

n  dz 


0  . 


(8) 


Substituting-  frca  equation  (7)  into  equation  (8)  yields  the 
governing  equation  or  vertical  equation. 


1-So-  +  X  2  N2(z)w*  =  0 
dz2  n 


(9a) 


subject  tc  boundary  conditions, 

w*  =  0  at  z  =  0  (Sea  Surface) 

w*  =  0  at  z  =  -H  (Sea  Bottom) 


(9b) 
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Equation  (6)  written  in  the  form. 


72(p„)t  -  +  0<P n>x  =  0 


(10a) 


represents  the  Rossby  wave  equation  for  which  solutions  of 
the  form 


pn  =  EXP[i(kx  +  iy  -  ot>]  ,  (10b) 

are  sought  where  k  and  Jl  are  wavenumbers  in  the  x  and  y 
direction,  respectively  and  <r  is  the  frequency.  Substitut¬ 
ing  from  (10b)  into  (10a)  yields  a  form  of  the  dispersion 
relationship  fcr  baroclinic  Rossby  waves,  viz.. 


-  8k  -  a (k*  + 


l2) 


of" 


(11) 


Equation  set  (9)  forms  a  Sturm-Liouville  eigenvalue 
problem  for  the  vertical  motion,  where  An  are  the  eigen¬ 
values  and  w*  the  corresponding  eigenfunctions *  A  solution 
is  sought  for  a  finite  number  of  inodes,  m.  These  modes  cor¬ 
respond  to  the  first  a  baroclinic  Rossby  wave  modes. 

If  S  (z)  were  independent  of  depth  at  a  fixed  loca¬ 
tion  and  time,  the  solution  for  wj  could  be  readily 
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functions  of 


determined  analytically  as  trigonometric 
(nirz)/H  (see  for  example  Eedlosky  (1979,  pp.  356-360)),  pro¬ 
vided 


n 


nir 

m 


n  •  192)3j«  • 


(12) 


In  certain  instances  an  analytical  solution  to  equation  set 
(9)  can  be  found  provided  tr\2)  is,  or  approximates,  a  well 
defined  function  (see  for  example  Sillmott  and  Mysak  ( 1980) 
who  solve  the  vertical  problem  with  N^z)  as  an  exponential 
function  cf  z)  .  Realistic  distributions  of  ii2are  depen¬ 
dent  upon  depth  and  generally  not  well  defined  functions  of 
z,  therefore  the  solution  of  equation  set  (9)  is  obtained 
numerically. 

2.  Determination  of  N2  and  The  Isotherm  Displacement 
Field  - 

The  Brunt-Vaisala  frequency  squared  is,  as  defined 
earlier,  equal  to  -(g/p0).pz  .  Since  the  model  neglects 

salinity,  the  simple  state  relation. 


P  =  P0U  -  a(T  -  T0)]  ,  (13) 

is  used  to  convert  the  Brunt-Vaisala  frequency  squared  from 
a  function  cf  the  vertical  mean  density  gradient  into  a 
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function  of  the  vertical  mean  temperature  gradient,  where 
T.  is  a  constant  reference  temperature  and  a  is  the  coeffi¬ 
cient  of  thermal  expansion.  Differentiating  equation  (13) 
with  respect  to  z  and  substituting  into  the  definition  of 
N2  gives, 

N2  =  og  T  .  (1**) 
z 

Equation  (14)  is  used  to  calculate  N2  from  the  mean 
monthly  temperature  fields  output  by  the  ocean  model  for  the 
last  ten  years,  1969-78,  of  the  model  run.  Only  the  temper¬ 
ature  field  in  the  region  18.3  -  56. 9N  and  159. IE  -  125W 
shown  as  the  enclosed  box  in  Figure  1  is  used  in  the  calcu¬ 
lation  of  N2  After  completion  of  the  numerical  differentia¬ 
tion  (right  side  of  equation  (14)),  a  ten-year  time  mean,  N2 
at  each  grid  point  and  at  19  vertical  levels  (mid-levels  of 
the  20  model  levels)  '  is  calculated.  Finally  N2is  averaged 
over  28  east-west  grid  points,  giving  a  N2  versus  depth  pro¬ 
file  for  each  latitude  (20  north-south  grid  points).  Justi¬ 
fication  for  the  east-west  averaging  is  based  upon  little 
variability  of  N2with  longitude  over  the  region  at  a  partic¬ 
ular  vertical  level.  The  variability  of  N2  with  latitude,  on 


the  ether  hand,  is  such  that  north-south  averaging  is  not 
justified.  Price  and  Nagaard  (  1980)  also  found  that  mean  N2 
profiles  varied  ouch  more  with  latitude  than  with  longitude, 
and  thus  the  numerical  model  data  are  consistent  with  the 
corresponding  data  taken  from  the  North  Pacific  Ocean.  Pro¬ 
files  of  (N2  jf^2  ,  averaged  over  longitude,  are  shown  in 
Figures  2  to  5  for  the  latitudes  of  the  grid  points  selected 
in  this  investigation  (triangular  symbols  on  Figure  1).  The 
solid  curves  cn  each  figure  are  drawn  to  the  values  of  (N2)'£ 
at  19  vertical  levels  determined  by  equation  (14).  The 
dashed  curves  are  drawn  to  values  of  (N2),,/2  calculated  at 
every  five  meters  from  z  =  0  to  z  =  -4000  m  using  interpola¬ 
tion.  The  interpolation  is  done  numerically  by  means  of 
Aitken's  scheme  of  Lagrange  interpolation  described  by 
Hildebrand  (1956,  pp.  49-50)  using  the  values  of  (N2)1/^2at  the 
19  vertical  levels  and  an  assigned  value  of  (N2  J1^2  at  -4000  m 
as  tabular  input.  At  latitude  48. 8N  and  54. 8N  the  assigned 
bottom  value  for  (ii2  was  given  the  identical  value  calcu¬ 
lated  at  level  19  (2830  m) .  At  the  other  selected  latitudes 
the  assigned  bottom  value  was  determined  by  linear  extrapo¬ 
lation  using  the  levels  below  1000  m  (levels  15  to  19). 


The  profiles  shew  a  pronounced  variability  with 
latitude.  Below  approximately  900  m,  (if2  Jr2  is  small  (less 
than  0.0C1  1/sec)  and  nearly  constant  (model  simulated 
isothermal  condition  in  the  deep  ocean).  At  latitudes  north 
of  47N  this  constancy  is  evident  at  depths  below  150  m.  The 
maximum  value  cf  (N2  )1/2  (maximum  vertical  temperature  gradi¬ 
ent)  occurs  at  level  2  (16.5  m)  at  latitudes  south  of  22N 
and  between  latitude  40N  and  50N.  At  the  ether  latitudes 
the  maximum  occurs  near  the  surface  (5  m).  A  secondary  max¬ 
imum  occurs  between  latitudes  30N  and  45N  at  a  depth  between 
200  m  and  400  ■  (inferred  location  of  the  mean  thermocline)  . 

The  second  field  tc  be  computed  using  the  model  sim¬ 
ulated  mean  monthly  temperature  fields  is  the  isotherm  dis¬ 
placement.  The  isotherm  displacement,  D(z,t),  is  the 
displacement  of  a  monthly  mean  isotherm  from  its  long-term 
average  level.  It  is  defined  mathematically  as  -T*/(T)Z 
where  T’  is  the  temperature  fluctuation,  T1  =  T  -  T.  T  and 
(T)z  are  the  long-term  (ten-year)  averaged  temperature  and 
vertical  temperature  gradient,  respectively.  In  the  numeri¬ 
cal  calculation  of  D(z,t)  the  assumption  is  made  that  the 
temperature  field  is  linear  in  z  between  model  levels.  The 
isotherm  displacement  versus  time  curves  at  19  vertical 
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Figure  4.  Profiles  of  (N2^at  Latitudes  42. 7N  and  48. 8N 
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levels  from  the  ten-year  model  simulation  are  shown  in 
Figures  6,7,  and  8  for  the  grid  point  at  30. 5N  139. iw. 

By  virtue  of  its  definition,  D(z,t)  positive,  in 
general,  implies  a  temperature  colder  than  the  long-term 
mean.  D(z,t)  negative  implies  a  temperature  warmer  than  the 
long-term  mean.  The  isotherm  displacement  curves,  as 
depicted  in  Figures  6  to  8,  reveal  a  pronounced  seasonal 
signal  (annual  period)  at  levels  above  71.5  m.  At  100  m 
(level  6)  and  below  the  seasonal  signal  becomes  less  well 
defined  and  losses  definition  with  depth.  The  seasonal  sig¬ 
nal  at  the  upper  levels  is  primarily  the  model  simulated 
mean  monthly  temperature  fields  response  to  zhe  seasonally 
varying  solar  radiation  provided  as  model  inpun.  Below  100 
m  (Figures  7  and  8)  a  low  frequency  (long  period)  signal  is 
evident.  The  isotherm  displacement  curves  at  the  other 
selected  grid  points  (not  shown)  ,  in  general,  show  similar 
characteristics  cf  an  upper  level  annual  cycle  and  low  fre¬ 
quency  signal  at  lower  levels.  The  significant  features, 
proceeding  from  south  to  north  at  longitude  155. 9W,  are 
that,  1)  the  depth  at  which  the  seasonal  signal  remains 
distinguishable  increases  reaching  1760  m  (level  17)  at 
54. 88,  2)  the  amplitude  cf  the  isotherm  displacement  at 
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mid-levels  increases  from  a  maximum  of  approximately  50  m  at 
18. 3N  to  5000  m  at  54. 8N.  This  increase  in  amplitude  is  a 
result  of  t.he  near  isothermal  conditions  evident  below  150  m 
at  the  northern  latitudes  (i.e.  (T)z  is  small  thereby 
increasing  the  value  of  D(z,t)). 

3.  The  Solution  cf  The  Vertical  Problem 

Equation  set  (9)  is  solved  numerically  using  an 
Interactive  Ordinary  Differential  Equation  (IODE)  solver 
which  uses  a  fourth-order  Bunge- Kutta  integration  technique. 
The  eigenvalues,  Xn,  are  found  by  means  of  the  "c-rcot 
method”  described  by  Preisendorf er  and  Gonzalez  ( 1  973)  .  In 
applying  the  ”c-rcot  method"  equation  (9a)  is  separated  into 
two  first-order  differential  equations. 


^S2n  =  4 
dz  v 

(15a) 

||  =  -  n2(z)»;  , 

(15b) 

A  trial  value  is  choosen  for  X,  and  equations  (15a)  and 
(15b)  are  integrated  using  IODE  starting  at  the  sea  bottom, 
z  =  -H  (4000  m) ,  with  boundary  values  of  wn'  =  0  and  <f>  =  1. 
As  the  integration  proceeds  the  zeros  of  the  function 
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function  c  (z)  =  w*(z)  -  $(z)Z(z')  are  noted.  In  the 

above  expression  Z(z')  equals  v'(z,)/0(z')  where  z'  is  the 
value  of  z  at  the  bottom  boundary  which  changes  as  the  inte¬ 
gration  continues.  Since  the  boundary  condition  states  that 
w*  (-H)  =  v'fz')  =  0,  then  Z(z)  must  equal  zero  and  the  c 

function  reduces  to  c  (z)  =  w*(z).  The  eigenvalue  problem 

thus  reduces  to  finding  values  of  z  where  w'=  0. 

For  example,  for  the  (N2)1/2  profile  at  30. 5N  and  a 
trial  value  of  X  =  2.5  sec/m,  c(z)  =  w*  (z)  =  0  when  z  = 

-2451m,  -liCSm,  -607m,  -326m,  and  -112m.  Varying  X  by  0.1 
sec/m,  noting  when  w'=  0,  and  tracing  out  the  points  results 
in  five  curves  (eigencur ves)  as  shown  in  Figure  9.  The 
first  five  eigenvalues  which  satisfy  the  boundary  condi¬ 
tions,  equations  (9b),  are  found  by  the  intersection  of  the 
eigencurves  with  the  z  =  0  line.  The  first  five  eigenvalues 
for  this  example  are  shown  to  the  right  of  the  ordinate  in 
Figure  9;  namely,  X,=  0.4735  sec/m,  X^s  0.8934  sec/m, 

X3  »  1.3380  sec/m,  X4  *  1.7814  sec/m,  Xg  =  2,)Q23  sec/m. 

Justification  for  limiting  the  calculation  cf  eigenvalues  to 
the  first  five  modes  will  be  made  apparent  in  the  next  chap¬ 
ter.  An  advantage  in  finding  the  eigenvalues  in  this  manner 
is  that  the  depth  to  the  sea  bottom  need  not  be  fixed.  If 
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the  depth  to  the  bottom  were,  for  example,  3  km  rather  than 
4  km,  the  eigenvalues  satisfying  the  same  given  boundary 
conditions  could  be  obtained  by  noting  the  value  of  at  the 
intersection  cf  the  eigencurves  and  the  z  =  -1000  m  line. 
The  reciprocal  eigenvalues,  cn  which  correspond  to  the 
internal  phase  speeds  in  the  absence  of  rotation  associated 
with  each  barcclinic  mode,  are  given  in  Table  I  for  the 
latitudes  of  each  selected  grid  point.  The  tabular  values 
of  cn  show  a  decrease  with  latitude  and  with  increasing  mode 
number. 

Using  expendable  bathythermograph  (XHT)  data  taken  at 
weather  station  November  (30n,  140R)  during  July  1966  to 
December  1969,  Emery  and  Hagaard  (1976)  obtained  cn  values 
of  c1  =  2.99  m/sec,  c2  =  1.41  m/sec,  c3  =  0.92  m/sec,  c4  = 
0.72  m/sec,  and  c5  =  0.51  m/sec.  These  values  are  slightly 
higher  than  the  values  given  in  Table  I  for  latitude  30. 5N 
(latitude  cf  licdel  grid  point  closest  to  the  latitude  cf 
weather  station  November)  .  The  difference  can  be  attributed 
to  the  deeper  depth  at  wearher  station  November  (4755  m)  and 
to  larger  values  for  N2(i.e.,  the  values  for  N2calculated 
from  the  model  simulated  temperature  field  is  smaller  than 
observed)  . 
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Eecipr 

LATITUDE 

TABLE  I 

coal  Eigenvalues  (m/s 

NODE 

1  2 - 

ec)  for 

NUMBER 

~  ~T 

Modes  1  to  5 

4 

5 

18.3N 

2.25  7  8 

1.2442 

0.8505 

0.6274 

0. 5049 

24 .4N 

2.2331 

1.1966 

0. 8202 

0.6099 

0.4922 

30. 5N 

2.11 19 

1.1193 

0.7474 

0. 5614 

0.4582 

36. 6N 

1.79  37 

C.9774 

0.6245 

0.4719 

0.3883 

42. 7N 

1.3289 

C.75  10 

0.4736 

0.3620 

0. 3077 

48.8N 

0.9299 

0.4831 

0.3294 

0.2629 

0.2264 

54. 8N 

0 . 55  1 1 

0.2598 

0.2153 

0. 1722 

0. 1356 

eigenfunct ions  are  then  normalized,  following  Emery  and 
Magaard  (1976) ,  as 


w  '  (z) 

rn(2)  s  — - - 2. -  ,  (16) 

[  J  N2(z)(wn' (z))2dz]1/2 
-H 


where  n  represents  the  mode  number.  Since  equation  set  (9) 
represents  a  Sturm-Licuville  system,  the  eigenfunctions  are 
orthogonal  with  respect  to  the  weight  function  N2  (z)  over 
the  vertical  column  z  *  -H  to  z  =0,  viz. , 

o 

/N2 ( z)w  ' (z)w  ' (z)dz  s  0  ;  n  *  m  (17) 

n  m 

-H 
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Therefore,  normalization  of  the  eigenfunctions  results  in 
the  crthonormal  39t,{yn}  ,  viz. . 

O 

-  f  j}2<z)Yn(2)Ym(z)  dz  =  5nm  *  (18) 

-H 

where  6,,,,,  is  the  Kronacker  delta  function.  A  more  detailed 
discussion  concerning  the  orthogonality  condition  for  the 
eigenfunctions  is  given  by  Willmott  and  Sysak  (1980) .  The 
normalized  eigenfunctions  for  the  selected  latitudes  are 
shown  in  Figures  10  through  16.  The  characteristics  of  the 
eigenfunctions  are  related  to  the  profile  of  N2  (z) .  At 
latitudes  south  of  40N  (Figures  10  to  13)  the  nodal  points 
of  the  normalized  eigenfunctions  are  concentrated  near  the 
upper  half  of  the  column  in  response  to  the  model  simulated 
upper  ocean  li2  maximum  (maximum  temperature  gradient).  At 
the  northern  latitudes  (Figures  15  and  16),  the  normalized 
eigerf unctions  are  very  nearly  sinusoidal,  particularly  at 
the  higher  mcdes.  This  corresponds  to  the  near  constant 
profile  of  ( N2  (Figures  4  and  5)  below  150  m. 

The  analysis  procedure  is  concluded  by  fitting  the 
normalized  eigenfunctions  to  the  simulated  isotherm  dis¬ 
placement,  D  (z, t) .  This  is  done  in  the  time  domain  by 
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Figure  11.  Normalized  Eigenfunctions  at  24.  4N 
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Figure  14.  Normalized  Eigenfunctions  at.  42. 7N 
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Figure  16.  Normalized  Eigenfunctions  at  54.  8N 
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making  use  cf  the  completeness  of  the  eigenfunctions  to 
equate  D(z,t)  to  the  sum  of  the  normalized  eigenfunctions 
multiplied  by  an  amplitude  coefficient  An,  viz. . 

m 

D(z,t)  =  ^  An(t)Yn(z)  '  <19> 

n=l 


The  coefficients,  An  (t)  ,  are  readily  obtained  by  making  use 
of  the  ort hcncrmali t y  condition.  Multiplying  each  side  cf 
equation  (19)  by  H27m and  then  integrating  from  z  =  -H  to  z  = 
0  results  in 

°  -  ° 

J  N2(z)Ym(z)D<z ,t)dz  s  ^  £An(t)  J  N2(z>Yn(z)Ym(z)dz]  .  (20) 

».  n=l 


Appl  icaticn  cf  the  othogorality  condition  (18)  to  equation 
(20)  yields 


An(t) 


O 


/ 


N2(z)Yn(z)D(z ,t)dz  . 


(21) 


Equation  (21)  is  used  to  compute  An(t)  for  the  first 
five  modes  for  each  selected  grid  point.  The  deviation  from 


the  medal  fit,  as  defined  by  Emery  and  Magaard  (1976)  ,eqn.(22). 


is  also  calculated.  The  deviation  term  quantitatively  shows 


how  well  the  model  simulated  isotherm  displacement  fluctua¬ 
tions  are  fit  by  the  first  five  baroclinic  modes.  A  value 
of  DIV(t)  *  0  implies  that  the  isotherm  displacement  fluctu¬ 
ations  can  be  fully  explained  by  the  first  five  baroclinic 
modes. 


w  5 

l  N2(z)[D(z,t)  -  ^  An(t)yn(z)]2  dz 


-H 


n=l 


f  N2(z)  D2(z,t)dz 


DEVCt) 


(22) 
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The  coefficients,  An(t),  and  the  deviation,  DEV(t)  ,  cal¬ 
culated  at  each  selected  grid  point  are  shown  in  Figures  17 
to  24.  At  lower  latitudes  {Figures  17  to  19  and  Figure  24) 
the  amplitude  cf  coefficient  A,  is  generally  larger  than  the 
coefficients  A2  tc  A5 indicating  the  importance  of  the  first 
mode  at  the  southern  latitudes.  At  latitude  24. 4N  (Figure 
18)  a  seasonal  signal  is  evident  for  coefficients  A4  and  A5. 
The  seasonal  signal  of  the  higher  modes  becomes  more  appar¬ 
ent  with  increasing  latitude  indicating  that  the  seasonally 
varying  model  input  parameters  appear  to  be  reflected  in  the 
higher  modes  at  higher  latitudes.  At  the  northern  most  grid 
point,  54. 6N  (Figure  23)  the  overall  amplitude  increases  an 
order  of  magnitude.  Additionally,  the  pattern  of  each  coef¬ 
ficient  appears  nearly  identical  with  only  slight  differ¬ 
ences  in  magnitude. 

The  deviation  term  at  latitudes  18. 3N  and  24. 4N  remains 
generally  telcw  50%  during  most  of  the  ten-year  model  run, 
although  scae  tine  periods  are  higher.  North  of  24. 4N  there 


is  a  general  increase  ever  the  ten-year  period 


in  the 
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Figure  18.  An  (t)  and  DEV  (t)  at  24. 4H  155. 9W 
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Figure  19.  An  (i)  and  DEV  (t)  at  30. 5H  155. 9W 
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Figur®  20.  An(t)  and  DBV  (t)  at  36. 6N  155. 9M 


MODE  NUMBER  3 


MODE  NUMBER  4 


MODE  NUMBER  5 


gure  23.  An(t)  and  DEV  (t)  at  54. 8N  155. 9W 
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deviation  curve  with  the  exception  of  the  two  northern  mcst 
grid  points  (Figures  22  and  23)  .  As  a  comparison,  the  first 
five  coefficients  and  the  deviation  term  computed  by  Emery 
and  Bagaard  (1976)  using  48  months  of  data  for  weather  sta¬ 
tion  November  are  shown  in  Figure  25.  The  results  shown  in 
Figures  24  (model  grid  point  closest  to  weather  station 
November)  and  25  compare  reasonably  well.  Both  show  the 
largest  amplitude  of  coefficaint  A1  .  The  48-month  mean 
deviation  computed  for  weather  station  November  was  32% 
while  the  ten-year  mean  deviation  is  38. 2%  for  model  grid 
point  10. 5N  139.  IN. 

To  examine  the  relative  contribution  made  by  each  baro- 
clinic  mode  tc  fitting  tte  simulated  isotherm  displacement 
fluctuaticns  (equation  ( 1 9) ) ,  the  ten-year  mean  deviation 
(equation  (22))  was  calculated  for  each  selected  grid  point 
adding  one  mode  at  a  time.  The  results  are  given  in  Table 
II.  Column  1  (under  modes)  gives  the  ten- year  mean 
deviation  when  only  the  first  baroclinic  mode  is  fit  to  the 
isotherm  displacement  fluctuations.  Column  2  gives  the 
ten-year  mean  deviation  when  modes  1  and  2  are  fitted  to  the 
isotherm  displacement  fluctuations  and  so  on.  The  tabular 
values  show  the  importance  of  mode  1  at  the  selected  grid 
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Figure  25.  An  (t)  ,  and  DE7  (t)  for  weather  station  November 
(from  Emery  ana  Sagaard,  1976) 

points  south  of  30. 5N  (i.e.,  over  40%  of  the  deviation  is 
accounted  for  ty  the  first  baroclinic  mode  alone)  .  The  con¬ 
tribution  made  ty  the  fifth  baroclinic  mode  in  reducing  the 
mean  deviation  is  less  than  6 %  at  the  three  southern  most 
grid  points,  and,  except  for  the  northern  most  grid  point, 
is  less  than  ICS  elsewhere.  Based  upon  these  results,  it  is 
apparent  that  adding  higher  baroclinic  modes  (modes  greater 
than  5)  would  add  little  to  reducing  the  deviation  term  com- 


TABLE 

II 

Dean 

Deviation 

(3)  at 

Selected  Gr 

id  Points 

GRID  FOIST 

MODES 

1 

1+2 

1  +  2+3 

1+2+3+4 

1+2+3+4+5 

18. 3N 

155. 9W 

56.3 

46.  0 

36.  4 

28.6 

26.8 

24.  4N 

155. 9W 

54.1 

45.  9 

39.0 

35.0 

31.6 

30. 5N 

155. 9W 

57.4 

54.2 

50.9 

46.  9 

41.0 

36. 6N 

155. 9W 

65.2 

55.7 

5  3.  3 

45.3 

38.8 

4  2.  7N 

155.9a 

96 . 2 

68.  1 

60.7 

54.  2 
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34.  5 
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26.6 
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6  3.6 
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The  mean  deviation  (column  5  of  Tabxe  II)  generally 
increases  with  latitude  up  to  42.7s  but  then  decreases  at 
the  two  highest  latitudes.  Additionally,  the  contribution 
of  mcdes  2  to  5  in  reducing  the  mean  deviation  is  signifi¬ 
cant  at  the  northern  most  points.  To  examine  the  relevance 
of  these  results,  the  cutoff  period  (minimum  period)  for  the 
propagation  of  free  Rossby  waves  of  mode  n  (see  Emery  and 
Sagaard  (1976))  given  by  T  =  (4irfAn)//3  was  evaluated.  The 
results  are  shewn  in  Figure  26.  The  dotted  line  delineates 
the  annual  period  and  the  dashed  line  the  five  year  period. 


PERIOD  (MONTHS) 


The  curves  cf  Figure  26  shew  the  following  results: 


•  Based  upon  the  theoretical  cutoff  oeriod  annual  oeriod 
baroclinic  free  Fossby  waves  may  exist  north  cf  2$N  but 
only  lode  1  can  occur. 


•  Baroclinic  free  Fossby  waves  at  the  annual  p  XO  d 

net  exist  north  of  approximately  36  S  (recall  Magaard 
(1983)  stated  that  first  mode  baroclinic  Fossby  waves 
were  primarily  evident  between  latitudes  30  to  40  North 
in  the  North  Pacific  Ccean)  . 

•  Baroclinic  free  Fossby  waves  at  the  five  year  period 
cannot  exist  north  cf  approximately  52N. 

The  results  given  above  are  significant.  Since  the  free 
Fossby  wave  analysis  was  conducted  using  only  ten  years  of 
model  simulated  temperature  data,  free  Fossby  waves  of  peri¬ 
ods  greater  than  approximately  five  years  can  not  be  ade¬ 
quately  resolved.  Therefore,  the  results  obtained  for  model 
grid  point  54.  £N  155. 9H  losses  significance.  Likewise,  at 
grid  point  48. 6N  155. 9W  the  theoretical  cutoff  period  curves 
show  that  only  first  mode  taroclinic  Fossby  waves  at  periods 
greater  than  approximately  three  years  can  exist.  There¬ 
fore,  the  low  deviation  at  latitude  48. 8N  resulting  prima¬ 
rily  from  modes  2  to  5  is  not  necessarily  due  to  taroclinic 
Fossby  waves.  The  opposite  may  be  said  of  the  selected  grid 
points  at  latitudes  19. 3N  to  36.  6N.  At  these  latitudes  tha 
first  baroclinic  mode  is  dominant  and  the  cutoff  period  is 
4.5  months  at  18. 3N  and  12.8  months  at  latitude  36. 6N.  At 
latitude  42.7N  the  cutoff  period  is  21.4  months.  However, 


only  the  first  and  second  modes  can  exist  at  periods  less 
than  five  years.  The  mean  deviation  is  high  when  only  the 
first  and  seccnd  modes  are  fitted  to  the  isotherm  displace¬ 
ment  fluctuations  (96.2%  and  88.1%  in  columns  1  and  2, 
respectively,  of  Table  II).  Thus  at  model  grid  point  42. 7N 
155.94  the  model  simulated  isotherm  displacement  fluctua¬ 
tions  cannot  be  explained  by  free  baroclinic  Rossby  waves  in 
the  time  domain. 

Thus  far  the  examination  of  baroclinic  free  Rossby  wave 
signature  in  the  model  simulated  temperature  field  (or  iso¬ 
therm  displacement)  has  been  limited  to  the  time  domain 
only.  In  order  to  conduct  an  analysis  of  baroclinic  free 
Rossby  wave  signature  in  the  space  domain  requires  carrying 
out  a  cross-sprsctral  fit  of  the  x  and  y  direction  wave  num¬ 
bers  at  the  appropriate  frequencies.  In  the  context  of  this 
study  a  cursory  examination  of  the  wave  number  vectors  will 
be  made.  The  locus  of  wave  number  vectors  for  mode  1  at 
appropriate  frequencies  for  latitudes  18. 3N  and  36. 6N  are 
shewn  in  Figures  27  and  28.  The  locus  of  wave  numbers  was 
determined  using  equation  (11)  for  three  different  frequen¬ 
cies. 
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Center  Circle  <r=  2.69  X.10 
Outer  Circle  <x  =  1.99  X  10 


-12.5-10.0  -7.5  -5.0^??S==^fe==5rS^5.O 


The  spacing  between  model  grid  points  utilized  in  this 
study  is  cr.  the  order  of  226  km  in  the  north/south  direction 
and  251  km  (at  36N)  in  the  east/west  direction.  The  shor¬ 
test  wavelength  resolvable  by  the  model  is  equal  to  twice 
the  spacing  between  grid  points,  which  in  this  case  is  on 
the  crder  cf  approximately  500  km.  Therefore,  Rossby  waves 
resolvable  by  the  model  must  have  wavelengths  greater  than 
or  egual  to  the  2nd  value,  where  ad  is  the  grid  spacing. 
The  wave  numbers  in  the  x  and  y  direction  which  correspond 
to  a  wavelength  which  is  greater  than  the  2nd  wavelength  are 
located  withir  the  dotted  rectangle  in  Figures  27  and  28. 
Stated  another  way,  wave  numbers  outside  of  the  boxea-in 
region  correspond  to  wavelengths  too  small  to  be  resolved  by 
the  medal. 

The  inner  circle  in  Figure  27  is  the  locus  of  wave  num¬ 
bers  found  when  a  frequency  corresponding  to  a  six-month 
period  is  used.  The  middle  circle  and  outer  crrcle  are  for 
frequencies  corresponding  to  periods  of  nine  months  and  12 
months,  respectively.  Note  that  due  to  the  different  scale 
of  each  axis  the  circles  appear  as  ellipses.  Portions  of 
all  three  circles  fall  within  the  enclosed  region  indicating 
that  wavelengths  greater  than  the  2dd  wavelength  can  occur 
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at  the  aforementioned  frequencies  (periods).  The  circles  in 
Pigure  28  are  for  frequencies  corresponding  to  periods  of  18 
nontbs  (inner  circle),  24  months  (middle  circle)  and  30 
months  (cuter  circle),  respectively.  It  is  noted  that  at 
36. 6S  only  waves  having  a  period  of  about  30  months  or 
greater  have  wavelengths  long  enough  to  be  resolvable  by  the 
model.  It  fcllcws  that  waves  of  4ad  or  8ad  would  have  even 
greater  periods.  Therefore,  increasing  latitude  requires  an 
increase  in  the  Scssby  wave  period  to  allow  for  wavelengths 
long  enough  to  be  resolved  by  the  model  grid  spacing.  Addi¬ 
tionally,  using  equation  (11),  it  may  be  shown  that  for 
wavelengths  that  are  resolved  by  the  model  grid,  k2+e24c  A2  f2 
and  thus  the  period  2  * /<x  must  be  larger  than  the  cutcff 
period.  At  high  latitudes  such  long  periods  would  not  be 
well  resolved  by  the  ten- year  simulation.  Thus,  even  though 
low  deviations  result  at  high  latitudes  (column  5  of  Table 
II)  the  isotherm  displacement  variability  cannot  be  attrib¬ 
uted  to  actual  Eossby  waves  with  any  degree  of  certainty. 


IV.  CONCLUSIONS  AND  HECOHMBN D AIION S 


Eased  upon  the  preceeding  analysis  the  following  results 


are  summarized: 


•  On  the  average,  free  baroclinic  Rossby  wave  signature 
in  the  vcdel  simulated  isotherm  displacement  fluctua¬ 
tions  is  evident  in  the  time  domain  for  the  selected 
grid  Feints  south  of  40N. 


•  At  the  selected  grid  points  south  of  40N  the  first 
baroclinic  node  shows  predominance. 


•  Eased  upon  the  locus  of  wave  number  ellipses,  Rossby 
waves  at  the  appropriate  periods  are  long  enough  to  be 
resolvable  within  the  model  grid  spacing. 


•  Eased  upon  the  deviation  term,  certain  times  within  the 
ten-year  lodel  run  period  indicate  that  the  first  five 
baroclinic  modes  do  not  adequately  fit  the  isotherm 
displacement  fluctuations  |i.e.,  the  deviation  is  high, 
therefore,  the  isotherm  displacement  fluctuations  can¬ 
not  be  fully  explained  by  baroclinic  free  Rossby 
wavas) . 


•  At  selected  grid  points  north  of  40N  the  analysis 
results  are  inconclusive. 


•  There  is  reasonably  qccd  agreement  between  the  analysis 
results  at  ucdel  grid  point  30. 5N  139. 1W  and  that  con¬ 
ducted  by  Emery  and  Hagaard  (1976)  for  weather  station 
November  (3 CN  14QW)  . 


The  results  stated  in  point  number  five  above,  are  a  product 


of  several  factors.  The  model  run  time  (ten  years)  sets  a 
constraint  upon  the  periods  that  may  be  adequately  resolved 
by  the  analysis  method.  The  cutoff  periods  of  the  selected 


grid  points  ncrth  of  40N  gradually  begin  to  approach  a 
five-year  period,  and  wavelengths  resolved  by  the  model  grid 
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result  in  periods  considerably  greater  than  the  cutoff  (min¬ 
imum)  period.  Periods  greater  than  five  years  car.  not  be 


adequately  resolved  by  tie  analysis  method.  A  longer  run 
time  is  required.  Secondly,  since  rhe  ocean  model  is  a 
"fresh”  water  model  the  effect  is  to  dimenish  the  magnitude 
of  the  Brunt  — vais&la  frequency,  particularly  at  higher 
latitudes.  This  results  in  eigenvalues  that  are  large  com¬ 
pared  to  "real"  ocean  values  ( A  is  inversely  proportional  to 
N)  .  The  larger  eigenvalues  then,  in  turn,  increase  the 
value  of  the  cutoff  period  more  than  what  occur  when  only 
the  Coriolis  parameter  is  considered.  Finally,  a  seasonal 
signal,  resulting  from  the  seasonally  varying  model  input 
parameters,  dominates  the  isotherm  displacement  fluctuations 
over  a  greater  portion  of  the  water  column  at  the  northern 
selected  grid  points.  Thus,  the  seasonal  signal,  which  is 
primarily  a  thermodynamic  effect,  contaminates  the  first 
five  modes. 

Eased  upon  these  results  several  recommendations  for 
further  research  may  be  set  forth: 

•  The  seasonal  signal  resulting  from  thermodynamic 
effects  should  be  removed.  This  could  be  achieved  by 
either  neglecting  the  model  layers  above  the  main  ther- 
■ccline  cr  by  applying  a  filter  technique  to  remove  the 
the  regular  seasonal  signal.  The  former  suggestion 
would  be  impractical  at  the  high  latitude  model  grid 
points,  however,  due  to  the  penetration  of  the  seasonal 
signal  dcwn  to  the  bottom  of  the  water  column. 


•  This  study  set  the  mean  current  equal  to  zero.  Further 
studies  snculd  include  the  mean  current  in  +he  analy¬ 
sis.  particularly  ip  the  region  of  the  California  cur¬ 
rent,  and  North  Pacific  Current.  A  study  by  Kang  and 
Magaard  (1979)  has  shown  that  the  North  Pacific  Current 
has  a  significant  impact  on  the  first  baroclinic  shear 
mode . 


•  Since  the  model  neglects  salinity  further  research 
should  study  the  feasability  of  incorporating  salinity. 
This  would  provide  fcr  more  realistic  density  profiles 
and  improve  the  analysis  results  at  tne  higher 
latitudes. 


•  The  mechanisms  which  produce  Rossby  waves  in  the  North 
Pacific  Ccean  is  still  not  fully  understood.  The  anal- 
ys.'s  has  shown  that  certain  time  periods  within  the 
moral  run  pericd  scuth  of  latitude  40 N  cannot  be 
explained  fully  by  the  first  five  baroclinic  modes, 
v  Patently  some  other  factor  is  involved.  Further 

research  could  be  done  to  isolate  any  peculiarities  in 
the  the  synoptic  wind  field  used  as  model  input  since 
the  six-nourly  synoptic  wind  was  the  only  varying 
parameter  from  one  year  to  the  next.  The  study  should 
also  he  extended  to  investigate  the  source  of  Rossby 
waves  (i.e.,  whether  they  result  from  large-amplitude 
variations  cf  the  wind  stress  curl  near  the  eastern 
boundary  as  hypothesized  by  White  and  Saur  (1981)  or 
whether  seme  other  mechanism  is  involved). 


•  All  further  studies  should  include  comparison  to  appro- 

Siate  observational  studies  when  available,  and 
erety  include  a  model  evaluation. 


This  study  has  teen  a  first  step  in  analyzing  the  low  fre¬ 


quency  dynamics  of  the  Haney  (1980)  ocean  model  below  the 
mixed  layer.  A  complete  understanding  of  the  dynamics  of 
the  large-scale  variability  of  the  ocean  temperature  field 


below  the  mixed  layer  (thermoci ine)  is  the  primary  goal. 


This  understanding  will  in  turn  lead  to  improvements  in  the 
ocean  model  and  its  diagnostic  and  prognostic  capabilities. 


the  ultimate  gcal 
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